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The results of calculations of the second interaction (cross) virial coefficient B (l w for water vapor 
and air, based on enhancement data obtained at NBS at 30, 40, and 50 °C, are presented. Comparisons 
are made with the results of calculations based on the enhancement data of Folitzer and Strebel, 
Webster, and Goff et al. and with the results of the theoretical calculations of Mason and Monchick 
and of Chaddock. An empirical equation is given for interpolation and extrapolation. 

The random (one standard deviation) uncertainty in the mean values of B aw, arising from the scatter 
of previously obtained NBS values of the enhancement factor, is estimated to range from 0.7 percent 
ai 30 °C to 1.4 percent of 50 °C. The estimated systematic uncertainties range from 4 percent at 30 °C 
to 6 percent at 50 °C, respectively. 
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1. Introduction 

The water vapor content of real gases, and of air in 
particular, saturated under known conditions of pres- 
sure and temperature, is not predicted adequately by 
ideal gas laws [l]. 1 In air, the deviations from ideality 
must be accounted for in order to obtain accuracies 
better than 0.5 percent at pressures as low as 900 
millibars [2]. On any isotherm, the saturated water 
vapor content increases with pressure to some, as yet, 
undetermined limit, in a manner consistent with that 
described by Haar and Sengers [3]. For example, at 
°C and 200 bars, the water vapor concentration or 
density in air is about twice that of the pure phase. 

This increase in water vapor content with pressure is 
the algebraic sum of the increases in apparent vapor 
pressure because of the superimposed pressure of air 
(the Poynting effect), and the Van der Waals type 
interactions between different molecular species, and 
the decrease in apparent vapor pressure due to the 
solution of the air in the liquid water (the Henry's law 
effect). Of the three the largest is due to the nonideality 
of the gas phase (the Van der Waals type interactions). 

Given a real gas equation of state of a water vapor- 
air mixture, say an equation expressed in virial form, 
it is possible to derive theoretically an expression 
for the saturation water vapor content of the gas 
mixture as a function of the mole fraction of the 
constituents, the parameters of state, and the virial 
coefficients. In such a formulation, the air-water 



1 Figures in brackets indicate the literature references at the end of this paper. 



interactions are in large part characterized by the 
second cross virial interaction term for the air and 
water molecules. Unfortunately there are no definitive 
values for the second cross interaction virial coeffi- 
cients, although previously obtained values [2] may 
be sufficient for some applications. 

Attempts have been made to derive second cross 
interaction virial coefficients from statistical me- 
chanics [4, 5]. However, the assumptions (form of 
potential) on which these derivations are based are 
far from exact so that the predicted coefficients are 
of limited value. What are needed are good experi- 
mental values of the second (cross) interaction coeffi- 
cients. Recently, Hyland and Wexler [6] at NBS have 
reported precise experimental values of enhancement 
factors for water vapor in C0 2 -free (CO2 content on 
the order of 2 ppm) air at temperatures of 30, 40, and 
50 °C. It is the purpose of this paper to use these new 
data, as well as the limited older data in the litera- 
ture, in order to compute values of the second inter- 
action virial coefficient. The formal derivations and 
basic experiments strictly apply to C0 2 -free air. This 
limitation is unimportant when using the B aw values 
in real air situations, as any errors introduced by the 
roughly 300 ppm of C0 2 molecules should be well 
within the limits caused by the uncertainties in B a w 

2. Theory 

2.1. General Considerations 

We will derive an equation which relates the inter- 
action virial coefficient to the enhancement factor, 
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thermodynamic parameters of state, and several 
physical constants. 2 

Initially, it should be stated that the treatment of 
air as a single-component gas leads to no theoretical 
inconsistencies, for the following reasons: (1) Over 
the experimental range of temperatures, we can ignore 
chemical reactions. (2) The experiments involved a 
continuous air flow over the liquid surface, so that 
the molar ratios of the air components are constant, 
for all experimental conditions of pressure and 
temperature, once equilibrium has been established. 
(3) We deal with only the chemical potential of the 
water in each phase, and since the required chemical 
potential difference in the liquid phase depends only 
on the mole fraction of water (see eq (19)) and not on 
the amount of dissolved species, it doesn't matter 
that the gas mixture dissolved in the water is of 
different composition than that of the gas-phase air. 
(This difference arises because the degree of absorp- 
tion in water varies from one air component to the 
next.) Thus, in this paper, it will be assumed that air 
acts as a single component substance with a known 
molecular weight. 

Let air be in thermodynamic equilibrium with a sur- 
face of the condensed water substance. The chemical 
potential in of each component in the gas phase is 
equal to that in the condensed phase. For our purposes 
we need only examine the chemical potential of the 
water, thus 



^ T ^KK) = ^ T > p KK) 



(l) 



where superscripts g and c designate the gaseous 
and condensed phases, T is the thermodynamic 
temperature, P the total system pressure, n 1 ^ (where 
k = g or c) is the number of moles of water, and n'; t 
is the number of moles of air. Water is a vapor in the 
gas phase, and either liquid or solid (ice) in the con- 
densed phase. 

In the equations which follow, the independent 
variables are always P, T 9 and n':\ whether stated 
explicitly or not. It will be understood that the vari- 
ables which do not explicitly appear are being held 
constant, so that subscripts are unnecessary for that 
purpose. 

Consider the difference in chemical potential in 
each phase between pressure states Pi and P>. We 
may write 



M&(ft)-/4(JM-MS,(l , i)--M&(iM. 



We now write 






(3) 



(see e.g., [7]) where G k is the total Gibbs free energy of 
either the gas or condensed phase. 

To obtain G k we use the thermodynamic relationship 



1 The theoretical basis for our work, which is reviewed below, is discussed in most 
good thermodynamics texts. It has been used in various derivations similar to ours, for 
example, by Haar and Sengers [3], and Goff and Bates [37]. 



yk = 



dP 



(4) 



where V k is the total volume of the phase. The Gibbs 
free energy difference between states P\ and P2, 
obtained by integrating eq (4), is 



G k (P 2 )-G k (P l ) 



-r 



V«dP. 



(5) 



Differentiating with respect to n 1 ^ it follows that 



071,., 



- ( r 



VHP 



Substituting eq (6) into (2) yields 



(6) 



(7) 



2.2. Gas Phase 

The equation of state of a gas may be expressed in 
virial form as a power series in reciprocal molar volume 



RT v v 2 



or as a power senes in pressure 



El 

RT 



= 1 + B'P + C'P 2 + 



(8a) 



(8b) 



where P is the total pressure, T is the absolute thermo- 
dynamic temperature, v is the molar volume, R is the 
gas constant, B and B' are second virial coefficients, 
and C and C" are third virial coefficients. In the en- 
hancement measurements considered in this paper, 
volume is not one of the experimental parameters. For 
this reason eq (8b) is used in our derivation. 

The virial coefficients of the pressure series are re- 
lated to those of the volume series by 



(2) and 



b- = -L 

RT 



C' = 



C-B* 

(RT) 2 



(8c) 
(8d) 



The virial coefficients are functions only of tempera- 
ture; those of eq (8a) are derivable from statistical 
mechanical relationships [8] if the form of the inter- 
molecular potential is known. (A very large "if", 
indeed! See, for example, Hanley and Klein [9]). The 
second virial coefficients may be considered to express 
the effects of interactions between two molecules, the 
third virial coefficients may be considered to express 
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the effects of interactions among three molecules, and 
so forth. The number of terms or coefficients necessary 
to adequately represent Pv/RT will depend on the gas 
(or gases) involved and the PvT state. 

If the gas under consideration is a mixture, then the 
coefficients B, C, B\ C\ etc., become mixture virial 
coefficients B m[x , C mix , etc., and can be written in 
terms of the mole fractions of the pure components, 
the virial coefficients of the pure components, and 
quantities called interaction (cross) virial coefficients. 
For a two-component gas mixture, and in particular for 
water vapor-air mixtures, statistical mechanics shows 
that [8] 

" mix X (tP (id ' "% aX w-t* aw ' X w D ww (-^a) 

^ mix X( ( Li (iaa ~r~ oXaX w L< aww ' oXqXw*^ aaw ' X W L> www V"*v 

where x u and x w are the mole fractions of air and water 
vapor, B (l a and B ww are the second virial coefficients 
for pure air and pure water vapor, and C (um and C www 
are the third virial coefficients for pure air and pure 
water vapor. B aw is the second interaction virial 
coefficient expressing the effects of interaction be- 
tween an air molecule and a water molecule. C (ia w is 
the third interaction virial coefficient expressing the 
effects of interaction between two air molecules and 
one water molecule whereas C a ww is the third inter- 
action virial coefficient expressing the effects of 
interaction between two water molecules and one air 
molecule. 

The mole fractions of air and water vapor are 
given by 



n lt 



X(i 



and 



n (l -t- n w 

n w 
n (l + n w 



(10a) 
(10b) 
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P_ 

RT 

(11) 



Now V mix =V !i ; therefore eq (11) may be substituted 
into the left-hand side of eq (7). We let Pi be the pure 
phase saturation vapor pressure e 8 (T) of the water 
substance and P-> be any other greater total pressure 
P. When the total pressure P reduces to e s , then 
x a = and x w =l. For the sake of simplicity, the 
superscript g will be deleted. After performing the 
integration and then the differentiation the left-hand 
side of eq (7) becomes 



r-r 

ldn w Je 



V ■ dP 

r mix u,i 



= RTln — - x*B aa P + 2B aw x-f 



+ {P- e, - xlP)B ww - (C aww - B% w ) -^ 



r (i + 2x n )(i-^ a ) 2 

t-'www ~ *^aaaX (l ~T DLiawwX a \l Xu) 



Zxj(\-2Xq) (l + 3* fl )(l-*q) 3 



3^4 
■ BaaBwwXl(l — 3x (l ) (1 — X„) + B 2 aa — - 



pi 

RT 



The molar volume ^ mix is related to the total volume 
^mix by 



*mix l^oT n w)Vmix 



(10c) 



+ \B a 



2x1(2 -3x„)P 2 
RT 



B u 



6X*(1 -Xa) 2 P* 



RT 



Bo 



where n a and n w are the number of moles of air and 
water vapor, respectively, in the total volume V mix . 
Substituting eqs (10), (9), (8c), and (8d) into (8b) we 
obtain 



*mix \ n « ' n w) D ' 



2n n n w 
n„ + n w 



B a 



Tin ~ 
ri 1 



B(, (l 



Tla + Uu 



B u 



(n a + n w )~ 

3n 2 a n w 

(n a + n w ) 2 



., ^aaa ■ 



3n a n 2 



(n a + n w ) 2 



,+ 



(n a + n w ) 2 



2x 2 l (l-Xq)(l-3xq)P 2 



XwP 

The ratio which appears in eq (12) will be called 

the "enhancement factor" and be designated by the 
symbol /. It has been variously called "the coefficient 
/," "the function /," and "the correction factor /" 
[2, 10, 11|. It is closely related to the "vapor con- 
centration enhancement" of Haar and Sengers [3]. 
Thus 



/= 



Xw£_ 



(13) 



The quantity x u P may be thought of as an "effective" 
vapor pressure of the water substance in a real gas 
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mixture analogous to a partial pressure in an ideal 
gas mixture. As P approaches e. s , x w approaches 
unity, and at P = e 8 and x w = 1, the enhancement factor 
/becomes unity. 



2.3. Condensed Phase 

We will confine our consideration of the condensed 
phase to the liquid state. Consider the right-hand side 
of eq (7) which expresses the difference in chemical 
potential for the water substance in the condensed 
phase between two pressure states P\ and Pi. As 
before, we let P\ be the pure phase saturation vapor 
pressure e s (T) of the water substance and P 2 be any 
other greater pressure P. Performing the differentia- 
tion yields 

/iS,(7\ P, »«, ng)-iL°,(T, e„ nj, 0) =^- \ V dP 

an wJe s 

(14) 

where V c is the total volume of the condensed phase 
(i.e., liquid water containing dissolved air), and where 
we have emphasized that at P\ = e s , n% = 0. Thus 



fi c Je s )=^ w (e s ) 



(15) 



where fxl v (e s ) is the chemical potential of the pure 
phase liquid water substance at pressure e s . 

Consider now the pure phase (single component) 
liquid water substance at the same pressure states P 
and e. s . The difference in chemical potential is 



V>M p )-p l w (es) 



Je s 



dP 



(16) 



where v x w is the molar volume of pure phase liquid 
water. 

Combining eqs (14), (15), and (16) yields 
fil,(P)-^ a (P)=^- f'VdP- \''v\ dP (17a) 



or 



-JLJ^ yc dP=f\l dP+^ w (P)-nl(P). (17b) 



The laws of dilute solutions [7] predict that for a solu- 
tion of dissolved air in water the chemical potential of 
the solvent, in this case liquid water, at pressure P 
differs from that of the pure phase of the same sub- 
stance under the same pressure by an amount given by 



lJL c w {P)-H l w {P)=RT\nx% + C 



where x c w is the mole fraction of water in the solution, 
and C is a corrective term expressing the excess 
chemical potential of a real solution over that predicted 
by the laws of ideal dilute solutions. 

Substituting eqs (18) and (17b) into (14) one obtains 
lx c w {P)-^ w {e s )= P vldP + RTlnx^+C. (19) 

Because x c w +x € a = l we may write 



lnx c w = \n (l-<). 



(20) 



For very dilute solutions 3 Henry's law may be used 
to calculate x c a , i.e., 



xc = k(T,P)x{>P 



(21) 



where k(T, P) is the Henry's law "constant" which is 
a function of T and P, and x° a is the mole fraction of 
air in the gas phase. 



The term C in eq (19) is given by [7| 

c=w{x< a y 



(22) 



where W is a function of the interaction energies of 
the molecular species in the solution. The term C 
can be shown [see appendix 1| to be negligible in its 
contribution to B (IW and so is dropped from further 
consideration. 



Equation (19) therefore becomes 
*£(**) -*£(*)= I" v\dP + RT\n (1-kxlP). 



(23) 



Kell and Whalley [12] have shown that the specific 
volume of the pure liquid phase of the water sub- 
stance can be represented with high accuracy by the 
following equation of state: 

V(t P) 5 3 

v, t : t 'r '+nvtf-w (24) 

K(£, latm) f=oPi 

where V(t,P) is the specific volume at temperature 
t (Celiuis) and pressure P, V(t, 1 atm) is the specific 
volume at temperature t and standard atmospheric 
pressure (1 atm), and P A is standard atmospheric 
pressure. Kell [13] has shown also that 



V(t, 1 atm) = 



l + bt 

5 

2 anf 



(25) 



(18) 



3 At pressures up to at least 200 bars and at temperatures from to 100 °C, the mole 
fraction of water in an equilibrium solution of dissolved air in water is very near unity. 
At 0° and 200 bars, x c = 0.997. 
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It therefore follows that 

vl=M».V(t,P) = 

Let 



i-0 j-l 



l + bt 



2 a «<" 



f 



i£dP-*(7\J»). 

Inserting eq (27) into (23) we obtain 

= K,(P)- M ;;,(e.,) = 



(26) 



(27) 



■ dn - 4 



r.nj 



g(r,P)+«rin(l-^«F) (28) 



2.4. Second Interaction Virial Coefficient 

By equating eq (12) and (28) and rearranging terms 
we obtain 



RT\nf=g(T, P)+RT\n {\-kx„P) 

+ BaaXJf - B W w (P - e s - X$P ) 



+ C, 



x^P 1 



-B] 



(KKI prp 

2RT 



+ C„ 



3x-(l-2x„)P 1 



3x-(l-Xa)P- 



' t> WWW 



D(UlBw 



— B'L 



2RT ""'"' RT 

(l + 2x„)(\-x (l )-P--ej 
2RT 

xt(l-3Xa)(l-X„)P s 

RT 

e?-(l+3*„)(l-* a ) 3 P 2 



+ B 

+ B W W 



au.[-2x|P- 



2RT 

2x 3 a (2-3x a )P* 



li. 



RT 



6xi(l-x„)-P- - 
RT 



B = C 



x jpi +c 

(1(1(1 T>rp ' Wl 



-C a 



3s;j(l-2s„)/* 
2RT 

3x, l ;(i— x„)P 1 

RT 
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Xa 
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D2 


e*-(l+3*«)(l- 


-Xa) 


>/>> 



ww 2 RT 

y=a + 6 
D = -2x i f 

2xl(2-3x„)P* 

O L> aa Dww 



6^(1-jc„) 2 P2 



RT 



<t> = D-8 

-2xl(l-x„) (\-3x„)P* 



€ = 



then 

so that 
B„ 



RT 



RTlnf=y+<t>B„w+eB? m 



(29c) 

(29d) 
(29e) 

(29f) 
(29g) 
(29h) 

(29i) 



2e 2 [ 



*)'■ 



(y-RT\nf)\ (30) 



1/2 



The quadratic term appearing in eq (29i) introduces 
a small correction into the calculated value of B a w 
Therefore the value of B aw obtained from eq (29i) 
without that term must be nearly the same as that cal- 
culated from eq (30). It may be shown that this condi- 
tion is satisfied only if the minus sign is used in front 
of the radical in eq (30). 

3. Sources of Data 
3.1. Virial Coefficients 



2xl{l-x a )(l-3x„)P> 

aw 



RT 



(29a) 



There are several sources of data for the second 
and third virial coefficients of air [14, 15, 16J. Probably 
the best and most up to date values oiB (m are those of 
For the sake of simplicity in nomenclature, all super- Sengers et al. [16] which cover the temperature range 
scripts have been omitted, it being understood, how- frorn 1°0 t0 1400 K. Sengers et al. assign a standard 
ever, that mole fractions x a and x w refer only to the gas deviation of 0.4 cm 3 /mol to B„ (l . We have chosen three 
phase. Let sigmas as our best estimate of the maximum syste- 

matic error, i.e., 1.2 cm 3 /mol. The Hilsenrath et al. 

/T D\ . DTI /i , D),D >« t 14 ^ ValUeS ° f Caaa are USed here * THeSe ran ^ fr ° m 

a = g(T, P)+RT\n (1 - kx„P) + B (l(l x;,P 90 to 1500 K. We have compared these values with 

those of Hall and Ibele [15] over temperatures from 
— B ww (r — e s — xf t P) (29b) to 100 °C, our range of interest. The two sets disagree 
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by about 10 percent. This disparity was arbitrarily 
increased to 15 percent and assigned as the estimated 
systematic uncertainty in C fma . 

The values given by Goff [2] are used for the second 
and third virial coefficients of water vapor. Wexler and 
Greenspan [17] recently have shown that from to 
100 °C the Goff values lead to highly precise correla- 
tions of theoretical and experimental values of satura- 
tion vapor pressure. Other second and third virial 
coefficients for water vapor are available but only those 
of Keyes [18] are both experimentally based and cover 
our experimental range of interest. Goff assigned toler- 
ances to his values, which he stated were equal to twice 
his estimated probable error. These were converted 
to one-sigma errors, and compared to the differences 
between corresponding values of Keyes and Goff. For 
B ww , the agreement between Keyes and Goff, below 
70 °C , is better than 2 percent, while the Goff one-sigma 
uncertainties range from 8 to 2 percent between 30 and 
60 °C. It was felt that the more conservative Goff uncer- 
tainties were appropriate here, and as estimates of the 
maximum systematic uncertainties in B ww we used 
three-sigma errors. In the case of C www below 70 °C, 
the Keyes and Goff values agree to about 45 percent, 
whereas the one-sigma errors computed from Goff 
range from 200 percent at 30 °C to 52 percent at 70 
°C. These are felt to be overly conservative. Thus the 
difference between the Keyes and Goff values was 
taken as an estimate of the systematic uncertainty 
in C wv)w» 

Mason and Monchick [5] and Hyland and Mason 
[19] give the only known values of C a <iw and C a ww<> 
respectively. 4 These are theoretical calculations based 
on the Lennard-Jones (12-6) potential, the first cover- 
ing the temperature range —80 to -h300 °C, the second 
from to 100 °C. For purposes of this analysis we 
ascribe an uncertainty of 50 percent in these values 
although there is no genuinely sound basis for this 
choice, and the errors may be larger. 

Interpolation equations for the various virial coeffi- 
cients are tabulated in appendix 2. 

3.2. Saturation Vapor Pressure 

The following equation, formulated by Wexler and 
Greenspan [17], is used to calculate e. s , the saturation 
vapor pressure of water: 



TABLE 1. Coefficients to vapor pressure formulation eq (31) a 



E 


- 7.51152 X10 3 


£, 


9.65389644X10' 


E> 


2.3998970 X 10"- 


E-, 


-1.1654551X10 5 


E 4 


- 1.2810336 X10- 8 


£5 


2.0998405 X 10-" 


B 


-1.2150799X10' 



a Units: e, s , pascals; T, kelvins. 

3.3. Function g(T, P) 

The coefficients a;,, a and b, given by Kell and 
Whalley [12] and Kell [13] are tabulated in table 2. 
These are needed in order to compute v x w using eq 
(26) which, in turn, permits the computation of g(T, P) 
using eq (27). Kell and Whalley estimated that the 
standard deviation ofjhe differences between observed 
values of V(T, P)IV(t, 1 atm) and eq (24) does not 
exceed 10 ppm. Kell estimated that V(t, 1 atm), 
given by eq (25) has a standard error of 10 ppm or less 
for temperatures up to 100 °C. 



}nes = ^EiT^ + B\nT 4 s 



(31) 





Table 2. Coefficients to eq (26) 


(i 




Coefficients ctij 


i \ 


1 


2 


3 



1 
2 
3 
4 
5 


-50.9769X10- 6 
3.71999 XIO 7 

- 7.01760 XIO" 9 
6.00227 x 10"" 

-3.09041 X10-' 3 
5.93416X10 ' 6 


8.2627 X10 ! » 
-1.3794X10-'° 

3.4032 XIO"' 2 
-3.6432 XIO' 4 

2.0836 X10' 6 
-4.1744 X10-' 9 


-9.109 X10-' 3 

2.626 XIO' 4 

-8.913 X10' 6 

11.467 XIO-' 8 
- 7.102 xlO- 20 

14.841 XIO 23 


n 


Coefficients a n 



1 
2 
3 
4 
5 


0.9998396 

1.8224944 X 10 2 
-7.922210 XIO- 6 
-5.544846 XIO 8 

1.497562X10-' 
-3.932952 XIO' 3 




Coefficient b 




1.8159725 X 10 " 2 



"Units: t, °C;P, bar; P A , 1.01325 bar; v\ cm 3 /mol. 



where T 4 h is the absolute temperature on the Interna- 
tional Practical Temperature Scale of 1948 [20, 21] 
and e s is expressed in pascals. 5 The coefficients Et 
and B are given in table 1. 



4 Incorrect values of C„ „•„• are listed in table 1 of the paper by Hyland and Mason. A 
B'^ w term was omitted in the calculation. Although the error propagates through the paper, 
the preferred C„ uu - values (calculated from viscosities) remain unaltered; the effect on the 
conclusions of the paper is negligible. 

5 Pascal =1 N/m- = 10- 5 bar = 10- mb = 7.50062 X 10 - mm Hg. 



3.4. Constants 

The solubility data of Winkler [22, 23] for air in 
water, as reported by Dorsey [24], was used to calculate 
the Henry's law "constant" k at standard atmospheric 
pressure. In the absence of any known air data on the 
pressure dependence of this "constant," it was assumed 
that the percentage changes in k between the same 
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Table 3. 


Henry s law "constants" k for nitrogen and air dissolved in water { 




Nitrogen ; ' 


Air 


Air e 


Temp. 


Pressure, bars 


Pressure, bars 


Pressure, bars 




1 


50 


100 


l b 


50 c 


100 d 


1 


50 


100 





18.65 


17.00 


15.62 


23.25 


21.19 


19.48 


23.20 


21.20 


19.49 


5 


16.42 


15.11 


13.96 


20.45 


18.82 


17.39 


20.53 


18.81 


17.37 


10 


14.89 


13.69 


12.70 


18.26 


16.80 


15.58 


18.30 


16.79 


15.58 


15 


13.51 


12.41 


11.60 


16.49 


15.15 


14.15 


16.46 


15.12 


14.10 


20 


12.52 


11.42 


10.72 


14.93 


13.62 


12.78 


14.96 


13.75 


12.88 


25 


11.61 


10.68 


10.01 


13.84 


12.73 


11.92 


13.76 


12.64 


11.90 


30 


10.88 


10.00 


9.46 


12.83 


11.79 


11.16 


12.79 


11.77 


11.13 


35 


10.23 


9.44 


8.97 


12.00 


10.08 


10.53 


12.02 


11.10 


10.52 


40 


9.67 


9.02 


8.57 


11.34 


10.58 


10.05 


11.40 


10.58 


10.05 


45 


9.27 


8.70 


8.25 


10.84 


10.17 


9.64 


10.87 


10.19 


9.68 


50 


8.90 


8.44 


8.00 


10.44 


9.90 


9.39 


10.39 


9.89 


9.38 



a Data based on solubilities in references [24, 25, 26, 27, 28, 29, 30]. 

h Data based on solubilities in references [22, 23, 24]. 

c Calculated by assuming that the percent difference for air is the same as for nitrogen 
between 50 and 1 bars. 

d Calculated by assuming that the percent difference for air is the same as for nitrogen 
between 100 and 1 bars. 

e Calculated from eq (32). 

f Units: A-, 10 ~ 6 mole fraction per bar; P, bars. 



pressure levels at given temperatures were the same 
for air as for nitrogen, for which there are solubility 
data both at atmospheric pressure [25, 26, 27, 28], and 
at high pressure for temperatures above 20 °C [24, 29, 
30]. Using solubility data for nitrogen in water and con- 
verting these to k, values for air were extropolated 
accordingly. Below 20 °C, solubilities of nitrogen in 
water were calculated using the method of Krichevsky 
and Kasarnovsky [31] and the fugacities reported by 
Demming and Shupe [32J. The polynomial 



10 6 A- = ]T dV 



(32) 



was fitted by the method of least squares to the Henry's 
law "constant" for air in water at 1, 50, and 100 bars, 
where / is in units of mole fraction per bar, C, are 
coefficients, and t is the temperature in degrees Celsius. 
Linear interpolations were employed between isobars. 
Values of k for nitrogen and air are given in table 3 
and for d in table 4. The values of k for air are con- 
sidered uncertain, at most, by 10 percent. 

On the unified carbon-12 scale the molecular weight 
of water M w is 18.0154 g/mol with a maximum total 



Table 4. Coefficients to eq (32) H 





Pressure, bars 




1 


50 


100 




1 

2 
3 


23.195 
-0.58037 

9.7392 xl 0- 3 
-6.5058 X10- 5 


21.197 
-0.51768 

8.2103 X10- 3 
-4.7585 X10- 5 


19.493 
-0.46024 

7.3781X10- 3 
-4.4401 x 10- 5 



"Units: t, degrees Celsius; k, mole fraction per bar. 



uncertainty of 0.0009 g/mol [33]. The gas constant R is 
equal to 83.1434 bar cm :j /mol K with a standard devia- 
tion of 0.0035 bar cm : 7mol K [34|. 

3.5. Enhancement Factor 

The NBS enhancement factor data [6] are given in 
table 5. It is estimated that the systematic uncertainty 
in /is 0.07 percent and the random uncertainty is 0.2 
percent. Two 30°C runs, reported as being suspicious 
in [6], have not been considered for the calculation 
of B (lw . 

4. Results 

Values of B (lw calculated from these data using 
eq (30) are given in table 5. The values were normalized 
from the experimental temperature to the nominal 
isotherm temperature. The changes, where they occur, 
are small. 

Our best estimates of the magnitudes of B (l w and 
associated random uncertainties are represented 
respectively by the means of the normalized isotherm 
values and the standard deviations of the mean nor- 
malized values. These are reported in table 6, along 
with the standard deviations of the individual 
determinations. 

5. Effect of Third Virial Terms 

Statistical mechanics predicts that the virial 
coefficients are functions solely of temperature [8]. 
An apparent dependence of B (m on pressure can be 
introduced via the calculations if a sufficient number of 
terms are not included in the truncated infinite 
series of the equation of state. If this equation were 
terminated at the second virial term rather than the 
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Table 5. Second interaction virial coefficients based 


on NBS enhancement data 






Saturation 


Saturation 


Total 






Interaction 


Normalized 




temp. 


vap. press. 


pressure 


Mole fractionl 
water vapor 


Enhancement 
factor 


virial coeff. 


inter, vir. coeff. 














Run 


°C 


mbar 


bar 






cm 3 /mol 


cm 3 /mol 


No. 
















t 


e s 


P 


x w 


/ 


Baw 


Baw 








Isotherm te 


mperature, 30 ° 


C 






39 


29.9994 


42.4278 


10.7312 


0.0040827 


1.0326 


-29.962 


- 29.962 


47 


30.0065 


42.4451 


15.1088 


.0029354 


1.0449 


-29.389 


-29.392 


42 


29.9999 


42.4290 


19.9344 


.0022534 


1.0587 


-29.115 


-29.115 


20 


30.0112 


42.4566 


30.1929 


.0015317 


1.0893 


-28.988 


- 28.944 


40 


29.9999 


42.4290 


35.5133 


.0013201 


1.1049 


-28.763 


-28.763 








Isotherm temperature, 40 c 


C 






34 


40.0104 


73.8143 


14.6417 


0.0052456 


1.0405 


-26.286 


-26.289 


19 


39.9969 


73.7612 


30.1957 


.0026497 


1.0847 


-26.978 


- 26.977 


27 


39.9648 


73.6350 


40.4557 


.0020231 


1.1115 


-26.124 


-26.114 


22 


40.1505 


74.3673 


50.3962 


.0016810 


1.1392 


-25.831 


-25.876 


21 


40.1366 


74.3122 


60.9945 


.0014285 


1.1725 


-26.158 


-26.199 








Isotherm temperature, 50 ° 


C 






32 


49.9898 


123.3233 


10.5512 


0.0120190 


1.0283 


-23.443 


- 23.440 


37 


49.9942 


123.3503 


10.6310 


.0119361 


1.0287 


-23.713 


-23.711 


36 


50.0042 


123.4115 


10.7441 


.0118214 


1.0292 


-23.915 


-23.916 


43 


50.0047 


123.4146 


10.7546 


.0117992 


1.0282 


-22.714 


-22.715 


46 


49.9971 


123.3680 


10.7953 


.0117446 


1.0277 


-21.977 


-21.976 


35 


50.0011 


123.3925 


10.7989 


.0117614 


1.0293 


-23.946 


-23.946 


31 


49.9891 


123.3191 


14.6231 


0.0087902 


1.0423 


-27.135 


-27.132 


33 


50.0064 


123.4250 


14.8620 


.0086263 


1.0387 


-23.399 


- 23.401 


30 


49.9967 


123.3656 


25.3585 


.0051773 


1.0642 


-23.186 


-23.185 


17 


50.0101 


123.4477 


30.1977 


.0044266 


1.0828 


-25.846 


-25.849 


16 


50.0050 


123.4164 


30.2324 


.0044135 


1.0811 


-25.117 


-25.118 


29 


50.0009 


123.3913 


36.5953 


.0036782 


1.0909 


-22.605 


-22.605 


26 


49.9766 


123.2426 


40.8462 


.0033228 


1.1013 


-22.502 


-22.495 


23 


50.0049 


123.4158 


50.3229 


.0027681 


1.1287 


-23.227 


-23.228 


24 


49.9784 


123.2536 


76.3139 


.0019376 


1.1997 


-23.088 


-23.082 


28 


49.9815 


123.2726 


96.9356 


.0016749 


1.2355 


-22.298 


-22.292 


25 


49.9788 


123.2561 


102.704 


.0015264 


1.2719 


-22.483 


-22.477 



Table 6. Mean normalized second interaction virial coefficients 
B aw based on NBS data 



Isotherm 

Temp. 


Mean 
normalized 

B aw 


Standard deviation 

of single 

determination 


Standard deviation 
of mean 


°C 


cm 3 /mol 


cm 3 /mol 


Percent 


cm 3 /mol 


Percent 


30 
40 
50 


-29.245 
-26.291 
-23.595 


0.46 

.41 

1.35 


1.6 
1.6 

5.7 


.21 
.18 
.33 


0.7 

.7 

1.4 



third, the terms associated with the parameters e, 
8, and /3 would reduce to zero and from eq (29i) we 
would have 

B aw =(RT\nf-a)ID. (33) 

The change introduced by using eq (33) instead of 
(30) is significant as shown graphically in figure 1 
using the NBS enhancement data. Because of this the 



third virial coefficients were included in the calcula- 
tion of B aw , even though there are large uncertainties 
associated with them. It will be shown that the accuracy 
in B (l w is relatively insensitive to the accuracies of the 
third virial coefficients. 

6. Error Analysis 

An analysis was made of the effect of suspected 
systematic and random errors on the accuracy of the 
values of Baw given in the results. 

6.1. Enhancement Factor 

It can be shown that the uncertainties in the experi- 
mental parameters P, e s , and x w contribute to the error 
in Baw primarily from their appearance in the enhance- 
ment factor term and negligibility because of their 
presence in the other terms of eq (30). 6 Therefore, for 



8 This was checked by making arbitrary changes in these parameters and noting the 
corresponding changes in the individual terms contributing to Baw, as well as in B au itself. 
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S -4.0 — 



50 60 70 

PRESSURE, BARS 



FIGURE 1. The percentage changes in calculated values o/B aw , when third virial 
terms are dropped from eq (30). 



purposes of error analysis, the higher order virial where x a is assumed roughly equal to unity. The 



coefficients may be negieeleu and eq (33) used with 
adequate accuracy for examining the effect of an error 
in / on Baw Differentiating B a w with respect to /, 
replacing the differential with finite increments, and 
considering only magnitudes, we obtain 



A£ aM 



RT RTAf 

~2Pxtf S 2Pf 



(34) 




40 60 80 

PRESSURE, BARS 



100 



FIGURE 2. Percentage errors in B aw as a function of temperature 
and pressure, arising from the estimated systematic error of 7 
parts in 10 4 in the enhancement factor. 



aysiciiiciiic uiiceiiainiy iii tuC iiu>^> vSxUcs Oi CHuancc 
ment factor is 0.07 percent [6]. By substituting this 
into eq (34), and using the mean values of B a w given 
in table 6 to convert AB a w to percent, the curves 
shown in figure 2 were generated. It may be noted that 
for a fixed relative uncertainty in /the corresponding 
uncertainty in B a w decreases with increasing pressure. 
This relationship suggests that for interaction virial 
coefficient determinations it would be advantageous 
to perform enhancement measurements at high 
pressures. However, as we will show below, this gain 
in accuracy is offset by the increasing uncertainties 
contributed by the terms in eq (30) containing the 
third virial coefficients. 



6.2. Virial Coefficients 

The effects of the estimated systematic uncertainties 
in the virial coefficients (given in sec. 3.1) on B (l w were 
calculated using eq (30). The results are given in table 
7 at the experimental pressures and temperatures. 



6.3. Henry's Law 

It is shown in appendix 1 that dropping the correc- 
tion term C to Henry's law, (eq (22)), may lead to an 
uncertainty in Baw of about 0.13 percent. An uncer- 
tainty is contributed also by the "constant" k. The 
second interaction virial coefficient may be written as 
follows: 

D -RTkxaP 

Kaw~ ~ p — + various terms. (35) 
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Table 7. Estimated systematic uncertainties in B a 





Press. 


Source of error 


Quadrature e 


Temp. 


Baa 


tf WW 


C« ana 


C WWW 


L> aaw 


Caww 


/ 


Other 


Virial 
















coefficient 


All 








Estimated error in parameter, percent 






errors 
only 


errors 




( a ) 


( b ) 


15 


( c ) 


50 


50 


.07 


( d ) 




°C 


bar 


Estimated error in B (IW , percent 


Estimated error in B aw , 
percent 


30 


10.7 


2.0 


1.9 


0.14 


0.05 


0.42 


0.42 


2.8 


0.14 


2.8 


4.0 




19.9 


2.0 


1.1 


.26 


.03 


.81 


.45 


1.5 


.14 


2.5 


2.9 




35.5 


2.0 


.68 


.47 


.02 


1.4 


.47 


0.85 


.14 


2.6 


2.7 


40 


14.6 


2.2 


1.6 


.20 


.08 


.62 


.66 


2.4 


.14 


2.9 


3.7 




40.5 


2.1 


.66 


.56 


.04 


1.7 


.71 


.85 


.14 


2.9 


3.0 




61.0 


2.1 


.48 


.84 


.03 


2.6 


.75 


.57 


.14 


3.6 


3.6 


50 


10.6 


2.6 


2.5 


.16 


.21 


0.46 


.98 


3.8 


.14 


3.8 


5.4 




25.4 


2.5 


1.1 


.38 


.10 


1.2 


1.0 


1.6 


.14 


3.2 


3.6 




50.3 


2.4 


0.64 


.76 


.06 


2.3 


1.1 


0.79 


.14 


3.6 


3.7 




76.3 


2.3 


.48 


1.2 


.04 


3.5 


1.2 


.52 


.14 


4.5 


4.6 




103 


2.3 


.40 


1.6 


.04 


4.8 


1.2 


.39 


.14 


5.7 


5.7 



a Estimated 3cr error is 1.2 cm 3 /mol; error is 17, 23, and 34 percent at 30, 40, and 50 °C, respectively. 
b Estimated error is 24, 16, and 11 percent at 30, 40, and 50 °C, respectively. 
c Estimated error is nominally 43 percent at 30, 40, and 50 °C. 

d Quadrature of errors contributed by Henry's law, gas constant, and ignoring the correction to the law of ideal 
solutions. 

e Computed by the square root of the sum of the squares of the individual terms. 



Hence, neglecting the sign, 



RT 
AB aw ~ — AA:- 13000 Mi. 



The estimated uncertainty in k is 10 percent. Over 
the pressure and temperature range of the NBS data, 
the variation in k is small. With sufficient accuracy for 
this calculation, A/c=l.lXl0 -6 and therefore AZ? au? 
= 0.014 cm 3 /mol. Assuming a nominal magnitude of 
26 cm 3 /mol for B aw , k induces an uncertainty in B aw 
of about 0.05 percent. 

6.4 Function g(P, 7); Gas Constant 

The largest contributions from the term g(P, T) to 
the uncertainties in B ( ,w arise from uncertainties in 
the specific volume of water V(t, P) and the molecular 
weight of water, while the largest contribution from the 
gas constant R arises through the enhancement factor 
term. The uncertainties in P and t in the NBS data are 
insignificant in their effect on V(t, P). The systematic 
errors in V(t, P), M w , and R (sees. 3.3 and 3.4) like- 
wise produce negligible uncertainties in the calculated 
values of B (lw . 

6.5 Estimated Systematic Uncertainty in B„ w 

The estimated systematic uncertainties in the 
individual parameters at the experimental values of 
P and T are summarized in table 7. These were com- 
bined by quadrature to give the estimated overall 



systematic uncertainty in B a w. Two quadrature 
columns are shown. The left column is the estimated 
uncertainty in the calculated value of B aw contributed 
(36) solely by the virial coefficient terms whereas the right 
column is the estimated uncertainty in B a w due to all 
suspected sources of systematic error. Along each 
isotherm, as the effect on B a w from the fixed percentage 
uncertainty in enhancement factor decreases with 
pressure, the effect from the uncertainties in the virial 
coefficients increases in such a way as to keep the 
overall systematic uncertainty in B aw more or less 
constant over the pressure range considered. The 
systematic uncertainty increases from 4 percent at 
30 °C to 6 percent at 50 °C, and applies to the mean 
value of Baw as well as to the individual values. 

6.6. Estimated Random Uncertainty in B aw 

The experimental standard deviations of the single 
determinations and of the mean value of B aw have 
been given (table 6) as our best estimates of the random 
uncertainties in B (lw . The random uncertainty in B (lw 
arises mainly from the random error in the enhance- 
ment factor. The latter, based on the residual standard 
deviations of fits of / to a pressure function, is 0.02 
percent at 30 °C, 0.13 percent at 40, and 0.26 percent 
at 50 °C [6] . The corresponding calculated uncertainty 
in a single B (lw determination, as a function of pressure 
along each isotherm, is indicated in table 8. Also given 
for each isotherm is a mean single-determination 
uncertainty, based on the three tabulated values. 
This calculated mean value should be comparable to 
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Table 8. Comparison of estimated random uncer- 
tainties in a single determination o/B aw 



Temp. 


P 


Random uncertainty in Baw 


Based on ran- 
dom error in f a 


Based on preci- 
cision of B llu 


°C 


Bar 


Percent 


Percent 


30 


10.7 
19.9 
35.5 


0.80 
.43 

.24 






Mean 0.49 


1.6 


40 


14.6 
40.5 
61.0 


4.4 

1.6 
1.1 






Mean 2.4 


1.6 


50 


10.6 
25.4 
50.3 
76.3 
103 


14.0 
5.8 
2.9 
1.9 
1.4 






Mean 5.2 


5.7 



a The standard deviation of a single determination 
of /is 0.02 percent at 30°C, .13 percent at 40°C, and 
0.26 percent at 50 °C. [6]. 



the experimentally determined precision of a single 
Baw determination. The latter values are repeated in 
table 8, and it is seen that the calculated and experi- 
mental random uncertainties are in reasonably good 
agreement. 

7. Comparisons 

There are three known experimental determinations 
of the enhancement of water vapor in air from which 



the interaction virial coefficient may be calculated. 
Politzer and Strebel [35) performed single saturation 
isotherm experiments at 50 and 70 °C at total pressure 
up to 200 bars. Webster [36] obtained values at— 35, 
—20, and 15 °C at total pressures up to 200 bars, 
also using the single isotherm saturation method. Goff 
et al. [37, 38, 39] measured a quantity closely related 
to the interaction virial coefficient, from 5 to 25 °C at 
total pressures near one bar. 

Values of B a w were computed from the Politzer and 
Strebel data using the same procedures, constants 
and virial coefficients that had been used with the 
NBS data. The 50 °C data of Politzer and Strebel 
yielded values of Baw that appear to scatter independ- 
ently of pressure around an average value, except for 
four points at 12 bars and below. These four points 
are suspect and so were discarded. The mean and the 
standard deviation of the mean are given in table 9. 
The 70 °C data of Politzer and Strebel show a strong 
monotonic pressure dependence, contrary to the 
predictions of theory. It is probable that there is a 
significant systematic error in the measurements 
although the source of this error is not obvious. Because 
of this, the 70 °C data were excluded from further 
consideration. 

Baw was computed similarly from Webster's data 
for each of his isotherms. At —35 and — 20 °C a minor 
change was introduced into the calculations to allow 
for the solid state of the condensed phase. This 
involved the use for the Henry's law "constant" the 
value for water at °C and 1 bar and for e s the appro- 
priate values of saturation vapor pressure with respect 
to ice [40]. The resultant error from the choice of k is 
negligible; in fact, within the uncertainty of the 
Webster measurements the Henry's law "constant" 
could be ignored. The mean value of B a w for each 
temperature, and the standard deviation of the mean 
are given in table 9. 



TABLE 9. Interaction virial coefficients from literature data 





Source 


Temperature 
°C 


Webster 


Politzer and Strebel 


Goff 


Mason and 
Monchick 


Chaddock b 


— Baw 

cm 3 /mol 


Baw 

cm'Vmol 


a 

- B„ w 
cm'Vmol 


B aw 

cm 3 /mol 


- Baw 

cm'Vmol 


°- R 

D aw 

cm : 7mol 


-Baw 

cm 3 /mol 


— Baw 

cm'Vmol 


-35 


64.91 
50.15 
38.01 


7.45 
3.96 
1.92 














-20 






48.2 
42.0 
39.3 
37.6 
36.8 
34.5 
32.3 
30.4 
28.5 


3.2 
3.0 
3.0 
2.9 
2.9 
2.8 
2.7 
2.7 
2.6 













50.4 
47.3 
45.4 
43.9 
41.0 
38.2 
35.8 
33.5 


62.0 


10 






58.1 


15 


34.91 


3.44 






56.1 


20 






54.4 


30 










51.1 


40 










48.1 


50 






23.16 


0.39 


45.2 


60 






42.5 















a Baw = mean value along an isotherm. 

b Neither Mason and Monchick nor Chaddock assigns an uncertainty to his theoretical values of B a w- 
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A set of Baw values were calculated by an empirical 
equation given by Goff [2], and converted to units 
consistent with those employed in this paper. The 
standard deviation ascribed to the Goff values were 
derived from the tolerances he assigned, which were 
stated to be two times the estimated probable error. 
The values are given in table 9. 

No attempt was made to estimate the systematic 
errors in the above values of B a w The standard devia- 
tions indicate only the precision of the measurements. 

Mason and Monchick [5] and Chaddock [4] have made 
statistical mechanical calculations of Baw by assuming 
that the forces between air and water vapor molecules 
can be represented by the Lennard-Jones (12-6) poten- 
tial and that the separate species potential parameters 
€o (the depth of the potential "well") and o"o (molecular 
"diameter") can be combined by simple combination 
rules. Their calculations have been extended to encom- 
pass our temperature range of interest and the values 
so obtained are given in table 9. 



experimental data [39]. The NBS values of B a w are 
smaller in absolute value than these other three sets of 
values. Although the NBS, Webster and the 50 °C 
Politzer and Strebel values appear to fall on a smooth 
curve, this may be fortuitous. Until there is additional 
experimental corroboration, this apparent consistency 
should be viewed with reservation. 

8. Smoothing Function 

Given experimental values of B aw over a reasonable 
temperature range, it is feasible to determine the form 
of the interaction potential between the water vapor 
and air species and therefore, to derive a theoretically 
based expression for interpolation and extrapolation 
of Baw Unfortunately the NBS data is too limited in 
temperature range to warrant this approach. 

A polynomial equation was fitted, therefore, to the 
NBS, Politzer and Strebel, and Webster data. The 
resultant expression is as follows: 
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FIGURE 3. B aw versus temperature, from various investigators. 



These several sets of values of B a w are compared 
graphically in figure 3. The curves representing the 
Goff, Mason and Monchick, and Chaddock values are 
similar in shape but displaced from each other. It 
is not surprising that the curves are alike for they are 
based on similar theoretical formulations, but with 
somewhat different potential constants. Even GofFs 
empirical equation is based on values calculated from a 
potential function similar to the Lennard-Jones that 
had been adjusted to yield results consistent with his 



B aw =^Diti 



(37) 



where D\ are empirical coefficients, t is the tempera- 
ture in deg Celsius and B aw is in units of cm 3 /mol. Di 
are listed in table 10, along with the residual standard 
deviation of the fit. The experimental and predicted 
values are compared in table 11. 

The differences between the experimental and pre- 
dicted values are well within the standard deviation of 
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the mean of the experimental values. The maximum 
deviation of this equation from NBS values of B aw is 
1.1 percent. 



Table 10. Coefficients to eq {37) 



D„ 


38.9221 


D, 


- .384587 


D 2 


+ .512266 x 10- 2 


D ;J 


-.117467 x 10-'* 


D 4 


+ .878093 x 10« 


0"res 


.842 cm 3 /mol 



where Z is a lattice coordination number (assuming a 
local crystal structure), No is Avogadro's number, 
€tj are interaction energies between the designated 
molecular species. 

Let the mole fractions x\ in eq (38) be converted to 
moles rii. Differentiating with respect to n w we obtain. 



dn w WX «' 



(40) 



This equation gives the excess chemical potential of 
the water with dissolved air over that predicted from 





Table 11. 


Comparison of experimental and predicted values o/B aw 




Temperature 


Interaction virial coefficient 


Difference 


Standard deviation 












of mean 




°c 


Experimental 


Predicted" 


Pred. — Exp. 


Experimental 


Source 


cm 3 /mol 


cm 3 /mol 


cm 3 /mol 


Percent 


cm 3 /mol 


Percent 




-35 


-64.91 


-65.01 


-0. 1 1 


-0.17 


7.45 


11.5 


Webster [36]. 


-20 


-50.15 


-49.74 


+0.41 


+0.82 


3.96 


7.9 


Webster ,36]. 





-38.01 


-38.92 


-0.91 


-2.34 


1.92 


5.1 


Webster [36]. 


15 


-34.91 


-33.95 


+0.96 


+ 2.75 


3.41 


9.8 


Webster [36]. 


20 




-32.48 












30 


-29.24 


-29.53 


-0.29 


-0.99 


0.21 


0.7 


NBS. 


40 


-26.29 


-26.47 


-0.17 


-0.65 


0.18 


0.7 


NBS. 


50 


-23.56 


-23.30 


+0.26 


+ 1.10 


0.33 


1,1 


NBS. 


50 


-23.16 


-23.30 


-0.14 


-0.60 


0.39 


1.7 


P&S[35J. 


60 




-20.30 













" Calculated using eq (37). 

Equation (37) is offered as a suitable smoothing 
function for Baw that is valid from 30 to 50 °C, the tem- 
perature range of the NBS experiments. We feel 
confident in using it from 10 to 60 °C. Because it fits 
the Webster data well, and within the uncertainty 
of that data, the equation may be used to —35 °C. How- 
ever, it should be noted that the systematic uncertainty 
in the Webster data is unknown so that the reliability 
of the predicted B a w in the temperature range below, 
say, 10 °C is unresolved. 

9. Appendix 1 

For a two component solution of equal-size mole- 
cules, such as air dissolved in water, it can be shown 
[7] that 

W{n a + n w )xaXw= AC^ (38) 

where AG^. x is the excess Gibb's function 7 , that is, the 

departure from that predicted through Raoult's law; 
n a and n w are the total number of moles of air and 
water; and x a and x w are mole fractions of air and 
water. The parameter Wis given by [7] 



Raoult's law and is therefore equivalent to the term 
C in eqs (18), (19), and (22). 

Mason and Monchick [5] and Chaddock [4] both give 
€aalk = 99.2 K and e W wlk = Sm K (where k is Boltz- 
mann's constant). Mason and Monchick give values of 
194 and 222 K for e aw lk, as determined from second 
virial coefficient data and viscosity data respectively 
whereas Chaddock gives 220.5 K. For the present 
purpose, assume e aw lk = 210 K. If a coordination num- 
ber of 15 (which should lead to a conservatively large 
estimate of C) is assumed, then W is —3.69 X 10 4 bar 
cm 3 /mol. Suppose that this is in error by a factor of 5, 
then W might be of the magnitude -1.84 X 10 5 bar 
cm 3 /mol. 

The mole fraction of air dissolved in water x a may 
be estimated from Henry's law using the constants in 
table 3. The worst case (i.e., the largest value of x a 
which, in turn maximizes C) occurs at 0°C and 100 
bars. Using this value for x a and ^=-1.84 X 10 5 atm 
cm 3 /mol, the magnitude of C is then 0.7 bar cm 3 /mol. 

Equation (30), with the C term included, can be 
approximated by 



r=^( 2€a „ 



i € ww ) 



(39) 



Baw ~ [Various terms] — 



g(T,P ) 
2P 



2P 



-™ (41) 



7 Lewis and Randall [7] (-all this the excess free energy of mixing and designate it by the 
symbol AF^. 



The lowest pressure present in the NBS enhancement 
data is 10 bars. The maximum contribution of C to B (lw , 
therefore, is about 0.035 cm 3 /mol. For B aw of the order 
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487-412 OL - 73 - 10 



of 26 cm 3 /mol, the relative contribution of C to B„ w is 
about 1.3 parts in 1000, which is about one order of 
magnitude smaller than the estimated uncertainty in 

D aw* 

10. Appendix 2 

The equations for the virial coefficients for water 
are essentially those of Goff [2]. B ww is given here with 
opposite sign to conform to eq (8a). C www is a modifica- 
tion of Goflfs equation to conform to eq (8a). Goff 
suggests the use of his equation for C www only at 
temperatures 60 to 100 °C whereas we extropolate it to 
— 20 °C. Although this is a potentially dangerous pro- 
cedure, no other alternative is available. 

The equations for the virial coefficients B (Ul and C (m a 
were obtained by fitting the Sengers et al. data [16] and 
the Hilsenrath et al. data [14] respectively, over the 
temperature range 260 to 380 K, to polynomials. 

The equations for the interaction virial coefficients 
Caaw and Caww were similarly obtained by fitting the 
Mason and Monchick [5] values from — 20 to + 120 °C 
and the Hyland and Mason [19] values from to 100 °C, 
respectively, to polynomials. 

The standard deviations of the fits for the latter four 
equations are given. These equations are; convenient 
empirical relationships and it is not intended that they 
be used beyond the fitted temperature ranges. 

rrw 72000 

B ww = 33.97 7f- X 10 ~r r cm : Vmol 



C www = 2.85558 -^ + B% w cm 6 /mol 2 

Baa = ~ 13.521 + 0.24234* - 0.10022 X 10 ^ 

+ 0.26880 X 10 -Vcm 3 /mol; a = 0.024 

Caaa = 1314.2 - 0.89453* - 0.36372 X 10 V 2 

+ - .48378 X 10-V - 0.18897 X 10V cm 6 /mol 2 ; 

o- = 0.17 

Caaw = 860.79 - 2.4203* + 0.92144 X 10" 2 * 2 

- 0.14568 X lO" 4 * 3 cm 6 /mol 2 ; a = 0.32 

C (t wwX lO 6 = -0.20263 + 0.52695 X 10 2 * 

-0.74761 xiOV 2 -f 0.57576 X10 6 * 3 
-0.18065 X lO 8 * 4 cm 8 /mol 2 ; a = 0.614 X 10- 3 

t= t + 273. 16, * = degrees Celsius. 
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